Role of surface roughness in superlubricity 

U. Tartaglino 

IFF, FZ-Jiilich, 52425 Julich, Germany and 

Democritos National Simulation Center, Via Beirut 2, 34014 Trieste, Italy 



o 
o 

(N 

u , 

Oh: 

< 



«4H 

o 



o 
o 



> 

en 

o 
o 



O 

o 



X 





(d) 




V.N. Samoilov 

IFF, FZ-Jiilich, 52425 Julich, Germany and 

Physics Faculty, Moscow State University, 117234 Moscow, Russia 

B.N.J. Persson 
IFF, FZ-Jiilich, 52425 Julich, Germany 

We study the sliding of elastic solids in adhesive contact with flat and rough interfaces. We 
consider the dependence of the sliding friction on the elastic modulus of the solids. For elastically 
hard solids with planar surfaces with incommensurate surface structures we observe extremely low 
friction (superlubricity), which very abruptly increases as the elastic modulus decreases. We show 
that even a relatively small surface roughness may completely kill the superlubricity state. 



I. INTRODUCTION 

Friction between solid surfaces is a common phe- 
nomenon in nature and of extreme importance in biology 
and technology T|. At the most fundamental level friction 
is (almost) always due to elastic instabilities at the slid- 
ing interface. At low sliding velocity an elastic instability 
first involves (slow) elastic loading, followed by a rapid 
rearrangement, where the speed of the rearrangement is 
much faster than, and unrelated to, the loading (or slid- 
ing) velocity. During the fast rearrangement the elastic 
energy gained during the loading phase is converted into 
irregular heat motion. The exact way of how the energy is 
"dissipated" has usually a negligible influence on the slid- 
ing friction force, assuming that the dissipation occurs so 
fast that no memory of it remains during the next elastic 
loading event. There are many possible origins of elastic 
instabilities, e.g., it may involve individual molecules or, 
more likely, group of molecules or "patches" at the inter- 
face which we have denoted by stress domains ,2^ 13,11, [1|. 
The most fundamental problem in sliding friction is to 
understand the physical origin and nature of the elastic 
instabilities. 

Elastic instabilities occur only if the lateral corruga- 
tion of the interaction potential between the solid walls 
is high enough, or the clastic modulus of the solids small 
enough. Roughly speaking, elastic instabilities can only 
occur if a characteristic elastic energy is smaller than a 
characteristic binding (or rather barrier height) energy. 
To understand this, consider the simple model illustrated 
in Fig. ^ In (a) a particle or atom is moving in a cor- 
rugated (substrate) potential. Connected to the particle 
there is a spring (spring constant k) which is pulled with 
the velocity v. If the spring is soft enough, or the po- 
tential barrier height U is high enough, i.e., U > ka? , 
the particle will perform stick-slip motion [Fig. QJa) , (b)] , 
involving slow elastic loading followed by rapid slip and 
dissipation of the (elastic) spring energy. In this case the 
(time averaged) force on the particle is independent of 




FIG. 1: (a) A particle moving in a corrugated potential. The 
particle is connected to a spring which is pulled with the ve- 
locity V. (b) Stick-slip motion with a soft spring, fca^ < U . 
(c) Continuous sliding with a stiff spring, ka^ > U . (d) Elas- 
tically stiff solid sliding on a rigid corrugated substrate, (e) 
Elastically soft solid sliding on a rigid corrugated substrate: 
here the atoms can rearrange to occupy binding positions. 



V. However, in the opposite case U < ka? [Fig. H^c)], 
the particle will follow the drive with a velocity which is 
always comparable to v. In particular, when the drive is 
on-top of the barrier so will the particle be [Fig.^c)]. In 
this case no rapid motion will occur and the (time aver- 
aged) friction force acting on the particle is proportional 
to V. 

In a more realistic situation one must consider the 
whole interface. In this case, depending on the elasticity 
and lateral barriers and the size of the contact area, elas- 
tic instabilities may or may not occur[y|. Assume first 
that an elastically very stiff solid slides on a rigid corru- 
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FIG. 2: An elastic bar with two binding sites (bump). The 
corresponding binding cavities of the substrate do not match 
exactly the positions of the binding sites, causing a compe- 
tition between the energy cost to stretch the elastic bar and 
the binding energy. 



gated substrate, [Fig.^d)]. In this case the atoms at the 
bottom surface cannot adjust to the corrugated substrate 
potential, and (for an incommensurate system) as some 
atoms move downhill other atoms move uphill in such a 
way that the total energy is constant. Thus, no elastic in- 
stabilities will occur during sliding, resulting in a very low 
sliding friction; this state has been termed superlubric 7] . 
However, when the block is elastically soft [Fig. H^e)], 
the atoms can rearrange themselves so that at any mo- 
ment in time almost all the atoms occupy positions close 
to the minima of the substrate potential. During sliding 
rapid jumps will occur from time to time where a particle 
changes potential well. In this case the friction is high 
and (at zero temperature) remains finite as the sliding 
velocity u — > 0. 

It is well known that elastically hard solids tend to 
exhibit smaller sliding friction than (elastically) soft 
materials |8|. One extreme example is diamond which un- 
der normal circumstances exhibits very low kinetic fric- 
tion coefficient, of the order of 0.01, when diamond is 
sliding on diamond. This can be explained by the nearly 
absence of elastic instabilities because of the elastic hard- 
ness of the material. However, if clean diamond is sliding 
on clean diamond in ultrahigh vacuum, a huge friction 
(friction coefficient of the order of /z « 10) is observed^. 
The reason is that the clean surfaces have dangling bonds 
(which are passivated by hydrogen and oxygen in the nor- 
mal atmosphere) so that the interaction between the two 
diamond surfaces is very strong and elastic instabilities 
(and wear processes) can occur resulting in a very large 
friction. 

It is important to note that even if solids are too stiff 
for elastic instabilities to occur on short length scale, the 
ratio between the effective elasticity and the amplitude of 
the lateral corrugation of the binding potential may de- 
crease when the system is studied at a longer length scale, 
which may make elastic instabilities possible on a longer 
length scale [lOl[ll| . To illustrate this, in Fig.|21we show 
a one-dimensional (ID) case, where an elastic bar (cross 
section area A) with two binding sites (bumps) is in con- 
tact with a substrate with two binding sites (cavities). 
When the binding sites on the elastic bar overlap with 



the binding sites (cavities) on the (rigid) substrate, the 
binding energy U is gained. In order to gain this binding 
energy the segment of the elastic bar between the bumps 
(length L) must elongate by the distance b. Thus the 
strain in the segment is b/L and the elastic energy stored 
in the elongated segment is Uc\ = VE{b/ LY /2 where the 
volume V = LA. Thus, f/d = AEb'^/{2L) which de- 
creases as the length of the segment L increases. It fol- 
lows that only when L > AEb^ / {2U) will the bound state 
have a lower energy than the non-bound state. Thus, 
only on a large enough length scale will the solid be elasti- 
cally soft enough for elastic instabilities to occur. In most 
practical cases one is not interested in a ID situation but 
rather in semi-infinite solids, which are intermediate be- 
tween the 2D and 3D case. For surfaces with randomly 
distributed binding centers this situation is much more 
complex than for the ID case because the effective elastic- 
ity changes as quickly with the lateral length scale as does 
the effective amplitude of the lateral corrugation of the 
binding potential (which from random walk arguments[l| 
scales as i)[l3,[lll. A detailed analysis of this situation 
indicates, however, that if no elastic instability can oc- 
cur at short length scale it is very unlikely that elastic 
instabilities will occur on any length scale of practical im- 
portance, except perhaps in the context of earthquakes 
[lOl [ll| . If instead of randomly distributed binding sites 
one assumes incommensurate surfaces, one would expect 
even weaker pinning effects, and it can be argued that 
in this case the ratio between the effective elasticity and 
the amplitude of the lateral corrugation of the binding 
potential increases as ^ L so that if no elastic instabili- 
ties occur at short length scale they cannot occur at any 
length scale [l3|- Below we will present numerical results 
where elastic instabilities do occur also for (nearly) in- 
commensurate structures, but in these cases one of the 
solids is elastically very soft so that instabilities can occur 
on a short length scale. 

The discussion above has focused on clean surfaces and 
zero temperature. Temperature is unlikely to have any 
drastic influence on superlubricity. However, it may have 
a strong influence on sliding dynamics when elastic insta- 
bilities occur. As soon as T > K, thermal noise is able 
to activate jumps over the barrier, i.e., to provoke pre- 
mature jumps before the (zero temperature) instability 
point is reached. It has been shown experimentally |a,[l3 
and theoretically|3, 13 that this has a crucial influence 
on friction dynamics at low sliding velocity. Similarly, 
weakly bound adsorbed atoms and molecules have a large 
influence on the sliding dynamics, and may strongly in- 
crease the friction force [I^ as the mobile adsorbates can 
adjust themselves in the corrugated potential between 
the block and the substrate, giving rise to strong pinning 
effects. In this paper we will not address the role of ad- 
sorbates or non-zero temperature, but we will focus on 
the simplest case of clean surfaces at zero temperature. 

Recently, superlubricity has been observed during slid- 
ing of graphite on graphite: in the experiment described 
in Ref. ^fl| a tungsten tip with a graphite flake attached 



to it is slid on an atomically flat graphite surface. When 
the flake is in registry with the substrate stick-shp mo- 
tion and large friction are observed. When the flake is 
rotated out of registry, the forces felt by the different 
atoms start to cancel each other out, causing the fric- 
tion force to nearly vanish, and the contact to become 
superlubric. 

Graphite and many other layered materials are excel- 
lent dry lubricants. The most likely reason for this is 
that the solid walls of the sliding objects get coated by 
graphite flakes or layers with different orientation so a 
large fraction of the graphite-graphite contacts will be in 
the superlubric state. This will lead to a strong reduc- 
tion in the average friction. However, the coated solid 
walls are unlikely to be perfectly flat and it is important 
to address how surface roughness may influence the su- 
perlubric state. In this paper we will show that even a 
relatively small surface roughness may kill the superlu- 
bric state. 

Lubrication by graphite flakes may even occur for 
diamond-like carbon (DLC) coatings, which may exhibit 
very low friction. Indeed Liu et al .17] have observed that 
a graphitized tribolayer is formed on top of diamond- 
like carbon coatings. Thus, also the excellent lubrication 
properties of DLC films might be caused by superlubric 
graphite contacts. We also note that DLC films are very 
hard and this will reduce the chance for elastic instabili- 
ties to occur jl^. 

In this paper we present atomistic Molecular Dynamics 
calculations of the sliding dynamics for contacting elas- 
tic solids with (nearly) incommensurate surface lattice 
structures. We consider both fiat and rough surfaces. 
We consider the dependence of the sliding friction on the 
elastic modulus of the solids. For elastically hard solids 
with flat surfaces and incommensurate surface structures 
we observe extremely low friction (superlubricity), which 
very abruptly increases as the elastic modulus is dimin- 
ished. We show that even a small surface roughness 
may completely kill the superlubric state. In order to 
study large systems we use a recently developed multi- 
scale approach 19j to contact mechanics where the num- 
ber of dynamical variables scales like ~ N"^ rather than 
as ~ N^ , where TV x iV is the number of atoms in the 
nominal contact area. 



II. MULTISCALE MOLECULAR DYNAMICS 

Let us discuss the minimum block size necessary in 
a computer simulation for an accurate description of the 
contact mechanics between two semi-infinite elastic solids 
with nominally fiat surfaces. Assume that the surface 
roughness power spectrum has a roll-off wavevector q — 
go corresponding to the roll-off wavelength Ao = 27r/qo- 
In this case the minimum block must extend Lx ~ Aq 
and Ly « Aq along the x and y directions. Furthermore, 
the block must extend at least a distance Lz ~ Aq in 
the direction perpendicular to the nominal contact area. 
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FIG. 3: Schematic structure of the model, (a) The fully 
atomistic model, (b) The multiscale smartblock model, where 
the solid in (a) is coarse grained by replacing groups of atoms 
with bigger "atoms". 



since the surface roughness with wavelength Aq affects the 
elastic block up to such a distance. Thus, the minimum 
block is a cube with the side L = Aq. 

As an example, if Ao corresponds to 1000 atomic spac- 
ings, one must at least consider a block with 1000 x 1000 
atoms within the xy contact plane, i.e., one would need 
to study the elastic deformations in a cubic block with 
at least 10^ atoms. However, it is possible to drasti- 
cally reduce the number of dynamical variables without 
loss of accuracy if one notes that an interfacial roughness 
with wavelength A will give rise to a deformation field 
in the block which extends a distance A into the solid, 
and which does not have any significant variation over 
distances much smaller than A. Thus when we study 
the deformation a distance z into the block we do not 
need to describe the solid on the atomistic level, but we 
can coarse-grain the solid by replacing groups of atoms 
with bigger "atoms" as indicated schematically in Fig. |31 
If there are N x N atoms in the nominal contact area 
one need n « IniV "atomic" layers in the z-direction. 
Moreover the number of atoms in each layer decreases in 
a geometric progression every time the coarse graining 
procedure is applied, so that the total number of par- 
ticles is of the order of N'^ instead of N^. This results 
in a huge reduction in the computational time for large 
systems. This multiscale approach may be implemented 
in various ways, and in Ref. J19l| we outline the proce- 
dure we have used in this paper which we refer to as the 
smartblock. 

The model presented above should accurately describe 
the deformations in the solids as long as the deformations 
vary slowly enough with time. However the phonons with 
short wavelength cannot propagate in the coarse grained 
region because the model does not implement enough de- 
grees of freedom for them. The short wavelength phonons 
get scattered back when they reach the coarse grained 
region, and this cannot be avoided within a standard 
molecular dynamics approach. However this is not a se- 
rious limitation for the present work: the static friction 



and the onset of sliding are still unaffected; moreover the 
superlubric sliding at low speed does not dissipate energy 
into phonons. On the other hand, if one had to study 
the spectrum of dissipated phonons, more advanced tech- 
niques to adsorb the energy of the phonons without back 
scattering have to be considered |2(1 l2l| . The only dissi- 
pation mechanism that we employed in our simulations is 
a Langevin thermostat at T = K (i.e., a viscous friction 
term) acting only on the atoms of the block far from the 
contact region. 

Figure O illustrates a case where the block is in the 
form of a cube with atomically flat surfaces. It is possi- 
ble to obtain curved surfaces of nearly arbitrary shape 
by "gluing" the upper surface of the block to a hard 
curved surface profile. This was described in detail in 
Ref. 1^. The elastic modulus and the shear modulus of 
the solid can be fixed at any value by proper choice of the 
elongation and bending spring constants for the springs 
connecting the atoms (see Refs. 22] and 19|). 

We note that with respect to contact mechanics, when 
the slopes of the surfaces are small, i.e. when the surfaces 
are almost horizontal, one of the two surfaces can be con- 
sidered flat, while the profile of the other surface has to be 
replaced by the difference of the two original profiles |2^. 
Thus, if the substrate has the profile z = /ii(x) and the 
block has the profile z = /i2(x), then we can replace 
the actual system with a fictive one where the block has 
an atomically smooth surface while the substrate profile 
/i(x) = /i2(x) — /ii(x). Furthermore, if the original solids 
have the elastic modulus E^ and i?2, and the Poisson ra- 
tio v\ and V2, then the substrate in the fictive system 
can be treated as rigid and the block as elastic with the 
elastic modulus E and Poisson ratio v chosen so that 
(1 - v^)lE = (1 - vDlE^ -f (1 - vDlE^. 

The results presented below have been obtained for an 
elastic flat block sliding on a rigid substrate. We consid- 
ered both flat and rough substrates. The atoms in the 
bottom layer of the block form a simple square lattice 
with lattice constant a. The lateral dimensions La: = ^xO- 



and Ly = Nya. For the block, N^, 



Ny = 48. Periodic 



boundary conditions are applied in the xy plane. The lat- 
eral size of the block is equal to that of the substrate, but 
for the latter we use different lattice constant b « a/(j>, 
where (j) — {1 + -\/5)/2 is the golden mean, in order to 
avoid the formation of commensurate structures at the 
interface. For the substrate, Nj; = Ny = 78. The mass 
of a block atom is 197 a.m.u. and the lattice spacing of 
the block is a = 2.6 A, so to get the same atomic mass 
and density of gold. The lattice spacing of the substrate 
is 6 = 1.6 A. We consider solid blocks with different 
Young's moduh from E = 0.2 GPa up to 1000 GPa. The 
Poisson ratio used for the block is ly — 0.3. 

The atoms at the interface between the block and the 
substrate interact with the potential 
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FIG. 4: The contact between an elastic block with a flat 
surface and a rough rigid substrate. Only the interfacial lay- 
ers of atoms are shown. The elastic modulus of the block is 
E — 100 GPa. The substrate is self-afline fractal with the 
root-mean-square roughness 3 A, fractal dimension Di — 2.2 
and roll-off wavevector qo = Sql, where ql ~ 2tt/Lx- The 
substrate and block interfacial atomic layers consisted of 
78 X 78 and 48 x 48 atoms, respectively. The applied pressure 
p = 10 GPa. Note the elastic deformation of the block, and 
that the real contact area is smaller than the nominal contact 
area. 



a — 1, Eq. (^ is the standard Lennard- Jones poten- 
tial. The parameter e is the binding energy between two 
atoms at the separation r — i^l^r^. When we study 
contact mechanics and friction without adhesion we put 
a = 0. In the calculations presented below we have used 
To = 3.28 A and e = 40 meV, which (when a = 1) 
gives an interfacial binding energy (per unit area) [22 

A7 « 4e/a2 « 23.7 meV/A^. 

As an illustration, in Fig. ^ we show the contact be- 
tween a flat elastic block (top) and a randomly rough 
rigid substrate (bottom) . Only the interfacial block and 
substrate atoms are shown. The substrate is self-affine 
fractal with the root-mean-square roughness 3 A (see 
Sec. lIII|) . Note the elastic deformation of the block, and 
that non-contact regions occur in the "deep" valleys of 
the substrate. Actually the real contact area is smaller 
than the nominal contact area. 



III. SELF-AFFINE FRACTAL SURFACES 

Consider a solid with a nominally flat surface. Let 
X, y,zhe a coordinate system with the x, y plane parallel 
to the surface plane. Assume that z = /i(x) describes the 
surface height profile, where x = [x, y) is the position 
vector within the surface plane. The most important 
property characterizing a randomly rough surface is the 
surface roughness power spectrum C(q) defined bvpa. 
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where r is the distance between a pair of atoms. When 



Here (...) stands for ensemble average and we have as- 
sumed that /i(x) is measured from the average surface 
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FIG. 5: Surface roughness power spectrum of a surface which 
is self-afiine fractal for qo < q < qi- The long-distance roll- 
off wavevector go and the short distance cut-off wavevector 
qi depend on the system under consideration. The slope of 
the logC — logg relation for q > qo determines the fractal 
exponent of the surface. The lateral size L of the surface (or 
of the studied surface region) determines the smallest possible 
wavevector qL ~ 2tt/L. 



plane so that (h) = Q. In what follows we will assume 
that the statistical properties of the surface are isotropic, 
in which case C{q) will only depend on the magnitude 
q = |q| of the wavevector q. 

Many surfaces tend to be nearly self-afhne fractal. A 
self-afRne fractal surface has the property that if part 
of the surface is magnified, with a magnification which 
in general is appropriately different in the perpendicu- 
lar direction to the surface as compared to the lateral 
directions, then the surface "looks the same", i.e., the 
statistical properties of the surface are invariant under 
this scale transformation|25J . For a self-affine surface the 
power spectrum has the power-law behavior 



C{q) ^ q 



-■2(H+1) 



where the Hurst exponent H is related to the fractal di- 
mension D[ of the surface via H = 3 — Df. Of course, 
for real surfaces this relation only holds in some finite 
wavevector region q^ < q < qi, and in a typical case 
C{q) has the form shown in Fig. |31 Note that in many 
cases there is a roll-off wavevector go below which C{q) 
is approximately constant. 

In our calculations we have used self-affine fractal sur- 
faces generated as outlined in Ref. [23. Thus, the surface 
height is written as 
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where, since /i(x) is real, i?(— q) = B{q) and (/>(— q) = 
— 0(q). If 0(q) are independent random variables, uni- 
formly distributed in the interval [0, 27r[, then one can 
easily show that higher order correlation functions in- 
volving /i(x) can be decomposed into a product of pair 
correlations, which implies that the height probability 
distribution P^ = {S{h — h{x))) is Gaussian ^5]. How- 
ever, such surfaces can have arbitrary surface roughness 



power spectrum. To prove this, substitute |j21l into (0) 
and use that 
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where Aq is the nominal surface area, gives 
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Thus, if we choose 

i?(q) = (2Vi)[C(q)]l/^ 



(4) 



where L — Aq , then the surface roughness profile (O 
has the surface roughness power spectrum C(q). If we 
assume that the statistical properties of the rough sur- 
face are isotropic, then C(q) — C{q) is a function of the 
magnitude q = |q|, but not of the direction of q. The 
randomly rough substrate surfaces used in our numeri- 
cal calculations were generated using © and (QJ and as- 
suming that the surface roughness power spectra have the 
form shown in Fig. [3 with the fractal dimension D{ = 2.2 
and the roll-off wavevector qo — 3qL, where qL — 2tt/Lx. 
We have chosen go — -^Ql rather than qo — q^ since the 
former value gives some self-averaging and less noisy nu- 
merical results. We also used qi = 27r/6 = 78go- The 
topography of the substrate with the root-mean-square 
roughness amplitude 3 A used in our numerical calcula- 
tions is shown in Fig. |B| 



IV. NUMERICAL RESULTS 

In this section we present the results of molecular dy- 
namics calculations of sliding of elastic blocks on rigid 
substrates. In all cases, unless otherwise stated, the up- 
per surface of the block moves with the velocity v = 
0.1 m/s, and the (nominal) squeezing pressure p is one 
tenth of the elastic modulus E of the block, i.e., p = OAE. 
The reason for choosing p proportional to E is twofold. 
First, we consider solids with elastic modulus which 
varies over several orders of magnitude, and it is not pos- 
sible to use a constant p as this would result in unphysical 
large variations in the elastic deformation of the block. 
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FIG. 6: Fractal surface of the substrate with the large cut-off 
wavevector qi = 2-K/b = 78 go- For a square 124.8 A x 124.8 A 
surface area. The fractal dimension D{ = 2.2 and the root- 
mean-square roughness amplitude is 3 A. 



FIG. 7: The static (curve 1) and kinetic (curve 2) friction 
coefficients as a function of the elastic modulus E of the block, 
for the flat substrate. In the calculation we have assumed 
the squeezing pressure cto — O.IE and the sliding velocity 
V — 0.1 m/s. 



Second, if two elastic solids are squeezed together with 
a given load, then as long as the area of real contact is 
small compared to the nominal contact area, the pres- 
sure in the contact areas will be proportional to the elas- 
tic modulus of the solids ^2|. Initially, when the block is 
pulled laterally it deforms loading clastic energy and the 
shear force between the surfaces increases gradually. The 
shear force reaches a maximum Fg at the onset of sliding. 
We used to such maximum shear force to calculate the 
static friction coefficient: ^s = Fs/Fn, -Fjv = pA being 
the normal force. During the sliding the shear force oscil- 
lates in time. We defined the kinetic friction coefficient 
as the ratio between the time averaged shear force and 
the normal load, i.e., fik = Fk/Fj^. 

Let us first assume that both the block and the sub- 
strate have atomically smooth surfaces. Fig.[7|shows the 
static and the kinetic friction coefficients as a function of 
the elastic modulus E of the block. Note the relatively 
abrupt decrease in the friction when the elastic modulus 
changes from Ei « 0.7 GPa to ^2 « 2 GPa. For E > E2 
practically no instabilities occur and the friction is ex- 
tremely small, while for E < Ei relatively strong elastic 
instabilities occur at the sliding interface, and the friction 
is high. For E = 0.2 GPa the static friction ^is > 2. This 
calculation illustrates that the transition from high fric- 
tion to superlubricity can be very abrupt; in the present 
case an increase in the elastic modulus by only a factor of 
~ 3 (from 0.7 to 2.1 GPa) decreases the kinetic friction 
by a factor of ~ 10^. 

In Fig.lHlwe show the time variation of the shear stress 
as a function of time when the elastic modulus of the 
block equals (a) E = 0.8 GPa and {h) E = 2 GPa. The 
elastic modulus of the stiffer solid is above the superlu- 
bricity threshold, and no (or negligible) elastic instabili- 
ties occur; the stress is a periodic function of time, with 
the period corresponding to the displacement 0.2 A. For 



the softer solid strong elastic instabilities occur during 
sliding, the shear stress is less regular (and the arrange- 
ment of the interfacial block atoms more disordered) than 
for the stiffer solid, and the (average) period is longer 
than 0.2 A. 

The regular pattern with period 0.2 A in Fig.jSfb) can 
be understood as follows. For our system, in the sliding 
direction there are 8 block atoms for every 13 substrate 
atoms. Assume first that the block (and the substrate) 
are perfectly stiff. In this case, the position of the 8 
block atoms will take 8 uniformly spaced positions within 
the substrate unit cell (lattice constant b), see Fig. 
Thus, a shift of the block with the distance 6/8 will take 
the system to a (geometrically) equivalent configuration. 
Hence, since 6 = 1.6 A we expect the periodicity of the 
shear stress to be b/8 = 0.2 A. When the block has 
a finite elasticity but above the superlubricity threshold, 
the atoms will relax somewhat in the substrate potential, 
but the configuration of the system will still repeat itself 
with the same period 6/8. 

Let us now consider the influence of surface roughness 
on the sliding dynamics. In Fig. 1101 we show the kinetic 
friction coefficients for an elastic block sliding on a rough 
substrate, as a function of the logarithm of elastic mod- 
ulus E of the block. The curves from top to bottom 
correspond to the substrate root-mean-square roughness 
amplitudes 3, 1, 0.3, 0.1 A and (flat substrate). For 
the substrate with the largest roughness, no superlubric- 
ity state can be observed for any clastic modulus up to 
E = 10^2 Pa. 

In Fig. ^] we show the kinetic friction coefficient as 
a function of the root-mean-square roughness amplitude 
of the substrate. The elastic modulus of the block 
E — 100 GPa. Note the strong decrease in the fric- 
tion when the root-mean-square roughness amplitude de- 
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FIG. 8: (a) The shear stress as a function of time for the fiat 
substrate. The elastic modulus of the blocli is E = 0.8 GPa. 
(b) The same as above but for the elastic modulus of the blocli 
E = 2 GPa. 
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FIG. 9: Commensurability ratio 8/13 between 1-dimensional 
chains. The upper image shows the period of 8 block atoms 
(light blue), i.e., 13 substrate atoms (yellow). The lower im- 
age is obtained by shifting the block for 1/13 of its lattice 
spacing. Block atoms occupy the same positions relatively to 
the substrate's hollows. 



creases below 0.3 A, which corresponds to a peak-to-peak 
roughness of roughly one atomic lattice spacing. 

Figure ITif a) shows the average (or nominal) shear 
stress as a function of time for the rough substrate with 
the root-mean-square roughness 3 A, and for the elastic 
modulus of the block E = 100, 50 and 20 GPa. Note that 
in addition to major slip events, several small slip events 
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FIG. 10; The kinetic friction coefficients for an elastic block 
sliding on rough substrates, as a function of the logarithm of 
elastic modulus E of the block. The curves 1-5 (from top to 
bottom) correspond to the root-mean-square roughness am- 
plitudes of the fractal substrate 3, 1, 0.3, 0.1 A and (flat 
substrate) . 
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FIG. 11: The kinetic friction coefficient as a function of the 
root-mean-square roughness amplitude of the substrate. The 
elastic modulus of the block E = 100 GPa. 



occur in all cases. These events correspond to local slip 
at some asperity contact regions before the major slip 
involving the whole contact area. In all cases, the time 
dependence of the shear stress remains periodic with the 
period 2.6 A, which corresponds to the lattice spacing of 
the block. Note also that for the elastically softer block 
{E — 20 GPa), the stress-noise increases after each ma- 
jor slip event; this is caused by the elastic waves (heat 
motion) excited during the (major) rapid slip events and 
not completely adsorbed by the thermostat. 

Fig. I12r b) shows the same as in (a) but now for the 
elastic modulus of the block E = 10 and 5 GPa. In this 
case the decrease of the elastic modulus of the block re- 
sults in the increase of both the static and kinetic friction. 
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FIG. 12: (a) The shear stress as a function of time for the 
rough substrate with root-mean-square roughness amplitude 
3 A. The elastic modulus of the block is E = 100, 50 and 20 
GPa. (b) The same as above but for the elastic modulus of 
the block E = 10 and 5 GPa. 



Fig. El shows the average displacement of the inter- 
facial atoms of the block (in the sliding direction) as a 
function of time. The root-mean-square roughness am- 
plitude for the substrate is 3 A. The elastic modulus 
of the block is E = 100, 10 and 5 GPa. Note that the 
slip distance for the major slip events increases as the 
elastic modulus of the block decreases, and that for the 
elastically hardest solid {E — 10^^ Pa) about a half of 
the forward displacement occurs between the major slip 
events. 

Fig.ll4lshows the average position of the interface block 
atoms in the z-direction (perpendicular to the sliding di- 
rection) as a function of time. Results are shown for 
the rough substrate with the root-mean-squarc rough- 
ness amplitude 3 A. The elastic modulus of the block 
is E = 100, 50, 20, 10 and 5 GPa. When the elastic 
modulus decreases, because of the adhesive interaction 
the block interfacial atoms come (on the average) more 
close to the rigid substrate, embracing the substrate as- 
perities. This increases the real area of contact between 
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FIG. 13: The average displacement of the interfacial atoms 
of the block as a function of time. The root-mean-square 
roughness amplitude for the substrate is 3 A. The elastic 
modulus of the block is _E = 100, 10 and 5 GPa. 
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FIG. 14: The average position of the interface block atoms in 
the 2-direction (perpendicular to the sliding direction), as a 
function of time during sliding. The root-mean-square rough- 
ness amplitude of the substrate is 3 A. The curves from top 
to bottom correspond to the elastic modulus of the block 
E = 100, 50, 20, 10 and 5 GPa. 



the surfaces and results in a higher friction. 

Fig. I15f a) shows the shear stress as a function of time 
for the rough substrate (root-mean-square amplitude 
3 A) and for the stiff block (£: = 100 GPa). The sohd 
curve is with the adhesion included, while the dashed 
curve is without the attractive part in the Lennard- Jones 
potential, i.e., with a = in Eq. ^. Note that with- 
out adhesion the major slip is not so pronounced as for 
the case with adhesion. Still the time dependence of the 
shear stress remains periodic with the same period 2.6 A, 
corresponding to the lattice spacing of the block. With- 
out adhesion, the shear stress curve is nearly symmetric 
around the zero-stress axis, and the kinetic friction coef- 
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FIG. 15: (a) The shear stress as a function of time calculated 
for the rough substrate (root-mean-square amphtude 3 A) 
and for the stiff block (E — 100 GPa). Curve 1 is calculated 
including adhesion, while Curve 2 is obtained without the 
attractive part in the Lennard- Jones potential, i.e., with a = 
in Eq. Q. (b) The average displacement of the interface 
block atoms (in the sliding direction) as a function of time for 
the same systems as in (a). 
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FIG. 16: (a) The shear stress as a function of time for the 
rough substrate (root-mean-square amplitude 3 A) and for 
the block elastic modulus _B = 10 GPa. Curve 1 is calculated 
including adhesion, while Curve 2 is obtained without the 
attractive part in the Lennard- Jones potential, i.e., with a = 
in Eq. Q. (b) The average displacement of the interface 
block atoms (in the sliding direction) as a function of time for 
the same systems as in (a). 



ficient (determined by the average shear stress divided by 
the squeezing pressure) is about 150 times smaller than 
when the adhesive interaction is included. 

In Fig. llSf b) we show the average displacement of the 
interface block atoms (in the sliding direction) as a func- 
tion of time for the same systems as in (a) . For the case 
without adhesion the major slip is not as abrupt as when 
adhesion is included. At every moment there is some 
lateral motion of the block interfacial atoms. 

Fig. I16f a) shows the shear stress as a function of 
time for the rough substrate (root-mean-square ampli- 
tude 3 A) and for the elastic block with the elastic mod- 
ulus E = W> GPa. The solid curve is with adhesion in- 
cluded, while the dashed curve is without the attractive 
part in the Lennard- Jones potential, i.e., with a = in 
Eq. Q. Fig. llBf b) shows the average displacement of the 
interface block atoms (in the sliding direction) as a func- 



tion of time for the same systems as in (a) . For the case 
without adhesion the major slip is not as abrupt as for 
the case with adhesion, and the sliding motion is nearly 
steady. In both cases, the time dependence of shear stress 
remains periodic with the period 2.6 A determined by the 
lattice spacing of the block. For the case with adhesion 
two small slips and a major slip can be observed in each 
period, and the kinetic friction is high. For the case with- 
out adhesion no elastic instability occurs, and the kinetic 
friction is very small. 

The roughness-induced increase of friction can be un- 
derstood by considering that the real contact involves 
only small regions, as shown in Fig. 01 The compensa- 
tion of the lateral forces that guarantees superlubricity 
between incommensurate surfaces (see Fig. ^ does not 
happen at the boundaries of the contact regions, neither 
it can happen on very small contacts with high curva- 



ture. Thus, the friction force should increase with in- 
creasing length of the boundaries between the contact 
and non-contact regions. This is completely different 
from what happens between commensurate walls or be- 
tween walls with very strong interactions, e.g., metals 
with cold-welded microjuctions, where the real area of 
contacts determines the friction force. 

The surfaces used in our simulations are self-afhne frac- 
tal up to the atomic scale, but in general the cut-off 
wavevector qi can be smaller, so that the typical size 
of the contact areas can be larger and the friction can 
be lower. In other words, the root-mean-square surface 
roughness alone is not enough to determine the amount 
of friction: a surface profile with wider mountains and 
valleys has to provide less friction than a surface pro- 
file dominated by short wavelength corrugation. In the 
same way, a surface with higher fractal dimension will 
have more roughness at the smallest wavelength, provid- 
ing higher friction for the same root-mean-square rough- 
ness and cut-off wavevector qi. In particular for high 
fractal dimensions the friction must depend dramatically 
on the cut-off wavevector. 



V. PRESSURE DEPENDENCE OF THE 
FRICTIONAL STRESS 

During sliding, the atoms at the sliding interface will 
experience energetic barriers derived from both the adhe- 
sive interaction between the atoms on the two opposing 
surfaces, and from the applied load. Thus, we may define 
an adhesion pressure Pad , and as long as Pad ^ P, where 
p is the pressure in the contact area derived from the 
external load, the frictional shear stress will be nearly 
independent of the applied load. Let us illustrate this 
with the system studied in Sec. IIVI Let us first consider 
the limiting case where the elastic modulus of the block is 
extremely small. In this case, in the initial pinned state 
(before sliding) all the block atoms will occupy hollow 
sites on the substrate, as indicated by atom A in Fig. 1171 
During sliding along the a;-direction, the atom A will 
move over the bridge position B and then "fall down" 
into the hollow position C (we assume overdamped mo- 
tion). The minimum energy for this process is given by 
the barrier height Se (the energy difference between the 
sites B and A) plus the work pa?5h against the exter- 
nal load, where a is the block lattice constant and 5h 
the change in the height between sites B and A (which 
depends on p). Thus the frictional shear stress Uf is de- 
termined by (Tf a^6 = (5e + pa^Sh, or 

cTf = Se/{ba'^)+pSh/b = {p^i+p)Sh/b, 

where we have defined the adhesion pressure pad ~ 
Se/{a^Sh). 

In our case (5e ss 3 meV and 6h « 0.008 A giving 
Pad ~ 10^° Pa. Thus, in the present case, only when the 
local pressure in the contact regions becomes of the order 
of ^ lOGPa, or more, it will start to influence the shear 




FIG. 17; A block atom moving (or jumping) from the hol- 
low site A over the bridge site B to the hollow site C. The 
maximum energy position along the trajectory is at site B. 



stress. This result is in accordance with our simulation 
results. Thus, for smooth surfaces, the shear stress acting 
on the block with the elastic modulus E = 0.5 GPa, 
squeezed against the substrate with the pressure p = 50 
and 150 MPa, is identical (~ 1 MPa) within the accuracy 
of the simulations. 

For inert materials such as rubber the adhesive pres- 
sure may be of similar magnitude as obtained above. 
Since the contact pressure for rubber in most cases is 
below 10 MPa, one may expect that the shear stress in 
the areas of real contact will be independent on the load. 
Recently, a strong dependence of the (apparent) shear 
stress on the squeezing pressure was observed for smooth 
Plexiglas balls sliding on very smooth silicon wafers cov- 
ered by silane lawrs|23. However, as one of us has ar- 
gued elsewhere [23, this does not reflect a fundamental 
dependence of the shear stress on the squeezing pressure, 
but has another origin. 



VI. SUMMARY AND CONCLUSION 

We have studied the sliding of elastic solids in adhe- 
sive contact with flat and rough interfaces. We consid- 
ered the dependence of the sliding friction on the elastic 
modulus of the solids. For elastically hard solids with 
planar surfaces with incommensurate surface structures 
we observe extremely low friction (super lubricity) , which 
very abruptly increases as the elastic modulus decreases. 
Thus, at the superlubricity threshold, an increase in the 
elastic modulus by a factor of ~ 3 resulted in the decrease 
in the frictional shear stress by a factor ^ 10^. We have 
shown that even a relatively small surface roughness may 
completely kill the superlubricity. For flat surfaces the 
shear stress is independent of the perpendicular (squeez- 
ing) pressure as long as the pressure p is below the adhe- 
sive pressure Pad, which typically is of the order of several 
GPa. 
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